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ABSTRACT 


We study the possibility that the universe has compact topologies T3, T? x R!, or S1 x R? using 
the seven-year WMAP data. The maximum likelihood 95% confidence intervals for the size L of the 
compact direction are 1.7 < L/Lo < 2.1, 1.8 < L/Lo < 2.0, 1.2 < L/Lo < 2.1 for the three cases, 
respectively, where Lo = 14.4 Gpc is the distance to the last scattering surface. An infinite universe is 
compatible with the data at 4.30. We find using a Bayesian analysis that the most probable universe 


has topology T? x Rt, with L/Lo = 1.9. 
Subject headings: 


1. INTRODUCTION 


General relativity is a local theory and does not pre- 
dict or constrain the global properties of the spacetime 
manifold describing our universe, which have to be con- 
strained through observations. Considerations of non- 
trivial spatial topology go back to as early as 1900 by 
Schwarzschild _(Schwartzschild (English transla- 
tion (Stewart et_alJ|1998)). In 1917 de Sitter noted that 
Einstein’s equations of general relativity did not con- 
strain the global structure of the spacetime 
[1917), while in 1924, Friedmann realized that non- 
positive curvature spaces could also have compact topol- 
ogy (Friedmannl (1924). A number of other papers were 


published in the following 


, and references therein) 
but the subject remained purely speculative un- 
til very recent accurate measurements of the CMB 


(cosmic microwave background) anisotropies by the 
COBE (Bennett_et_al||1996)) and WMAP (Jarosik et al. 
2011 2011 2011) satellites, 


and other observational advances in astrophysics and cos- 
mology. 

Various models for the topology of the universe (for 
the classification of different possible topologies see, 
e.g. (Barrow & Kodama [2001)) have been extensively 
studied recently and compared to the experimental data. 
The two most important ones are the Poincaré dodeca- 
hedral space and the 3-torus T?; these models are in best 
agreement with the experimental data. The Poincaré do- 
decahedral space arises by slicing the 3-sphere $° and 
thus has positive curvature, while the 3-torus is ob- 
tained by slicing infinite Euclidean space R? and there- 
fore is flat. Theoretical arguments about quantum cre- 
ation of the universe favor the flat case. Based on 
the Wheeler-DeWitt equation, Linde has argued 
(see also (Zeldovich & Starobinsky; i983) that 
compact flat universes are much easier to create than 
other models, and can naturally provide initial conditions 
for the onset of inflation. Furthermore, Roukema con- 
structed a measure on the set of compact manifolds and 
showed that non-flat models almost never occur while 


flat models occur almost certainly (Roukema & Blanloei 


(see also for a discussion of 
the Poincaré dodecahedral space versus the 3-torus). 
Following these arguments, we will focus our atten- 
tion entirely on the flat case in this work. Details 
on the Poincaré dodecahedral space and experimental 
data analysis are given in (Luminet et_al.[2003; 
2009; 2004) |2008; [Caillerie et al. 
Lew & Roukemal |2008; [Aurich et_al|[2005b), and other 
topologies in (Aurich & Lustig| |2010; [Niarchou & Jaffe 
2007; Amih et al peod Locat Creswell et ah 
Aurich et al.|(2005ġ [Roukemal 2005; 2004; 
Riazuelo et al] |2004 [Gomero et al- 
2000 2005 ! 
We analyze three different flat topologies in this work, 
My = T?, Ma = T? x R!, and ™ = St x R?, where 
the subscript denotes the number of non-compact direc- 
tions. We will generically refer to all three cases as a 
torus. Usual flat space is œ% = R. The topology 
of the 3-torus T? can be obtained by identifying the 
opposite edges of a parallelepiped. We only consider 
the simplest case of a rectangular parallelepiped with 
equal side lengths L. This has the highest number of 
symmetries which helps reduce the computational time. 
Moreover, it has been argued in that 
only in well-proportioned spaces is the quadrupole of the 
CMB temperature-temperature correlation function sup- 
pressed compared to the infinite universe. The surpris- 
ingly low observed quadrupole is one of the motivations 
to invoke a compact topology. For our case of the 3- 
torus, well-proportioned means that all three sides should 
be approximately equal. The topologies of the spaces 
T? x R! and S! x R? are obtained by compactifying only 
2 or 1 dimension respectively. Again, following the argu- 
ment of (Weeks et_al./2004)) and for the sake of simplicity, 
we consider only the case where the compactified dimen- 
sions of T? x R! have the same size. The size of the 
compact directions will be denoted by L. As L —> oo all 
three manifolds reduce to infinite flat space R3. 
Different approaches have been proposed for extract- 
ing information about the topology of the universe from 
the experimental data, the two most important ones be- 
ing the circles-in-the-sky test and the analysis of the 
CMB power spectrum. The basic idea of the circles-in- 
the sky test is that if the global structure of the space 
is smaller than the distance to the LSS (last scattering 


2 


surface), then the LSS will self-intersect in circles, pro- 
ducing correlations between circles with different centers. 
The detection of such circles can reveal the global prop- 
erties of space (for detailed description of the method 
ee omer 
2004alb)). The main disadvantage of the method 
is that it cannot be used if the size of the universe is 
bigger than the observable part of it (the distance to 
LSS). The one-year WMAP data has been analyzed with 
this method for signatures of non-trivial spatial topol- 
ogy (Cauihet all@00d), aae oub the possilility of 
compact spaces with a length-scale smaller than 24 Gpc. 
This limit has been extended by about 10% by the au- 
thors of 2007) who have also ruled 
out the possibility of Poincaré dodecahedral space. The 
authors of (Bielewicz & Banday|2011) have analyzed the 
most recent seven-year WMAP data with this method, 
putting a lower bound of about 27.9 Gpc on the size of 
the fundamental domain for a flat universe. 

The low-/ (i.e. large scale) portion of CMB correlations 
is sensitive to the topology of space, which gives rise to 
another method for detecting the topology. The torus 
preserves the homogeneity of infinite space but breaks 
rotational invariance. This implies that the power spec- 
trum of CMB temperature-temperature correlations does 
not contain all the possible information since the off- 
diagonal elements of the covariance matrix in the spher- 
ical harmonics expansion are non-zero in general, while 
the diagonal elements with equal / and different m values 
are not all equal to each other (see section |3|for more de- 
tails). Moreover, in it has been 
argued that the off-diagonal elements contain more infor- 
mation than the diagonal ones if the side length of the 
torus is less than twice the distance to the last scattering 
surface. Therefore, to gain all the possible information 
from the correlations of CMB anisotropies, one has to 
consider the full covariance matrix rather than just the 
power spectrum. 

The CMB correlation functions have been previousl 
used to analyze COBE (de Oliviera Costa & Smoo 
WMAP 
and three-year WMAP 2008 
The lower bound on the side length L of T? ob- 
tained from COBE data 
is L > 4.32h-'Gpc at 95% confidence, and 
L > 5.88h7-!Gpc at 68% confidence. For T? x R! 
and St x R? the lower bound obtained from 
COBE is L > 3.0h—! Gpe 
at 95% confidence. The authors of (Phillips & Kogut 
have obtained higher bounds for T?; L > 1.2 Lo 
at 95% confidence and L > 2.1 Lo at 68% confidence.! 
They have also found that the maximum likelihood 
occurs for L = 2.1 Lo (29Gpc). Several models of 
tori with different side lengths have been considered 
in with the conclusions L > 19.3 Gpe 
for T? and L > 14.4Gpe for S! x R?. The main result 
of is that the 3-torus with volume 
=~ 5 x 10% Gpc? (which corresponds to side length of 
17 Gpc) is well-compatible with the WMAP three-year 


1 We will give lengths in terms of Lo = 14.4Gpc, the dis- 
tance to the last scattering surface. The Hubble length is Hy = 
2.998/ho = 4.266 (0.703/ho) Gpc. 


data. The WMAP seven-year data has been analyzed 
for detecting signatures of the so-called half-turn space 
(lie aly “difeencs of Vie 
half-turn space from the 3-torus is that one of the edges 
is turned by 180° before identifying with the opposite 
edge) where the case of the 3-torus is also considered. 
Out of these works, only in 
and has the full covariance matrix 
been analyzed. For some earlier results on these 
topologies see also 
and references therein. 

Other methods for experimental detection of non- 
trivial topologies have been proposed and used to an- 
alyze the experimental data. The so-called spatial cross- 
correlation function method has been used in (Aurich 
to analyze the five-year WMAP data for signa- 
tures of a 3-torus topology with a result L = 3.85 Hg. 
They have also specified the orientation of the torus. 
In (Bielewicz & Riazuelo multipole vectors have 
been used to analyze the five-year WMAP data with 
the conclusion that a torus topology is slightly preferred. 
The authors of have looked 
for evidence for a 3-torus topology by the alignment of 
distant objects. They have put a lower bound on the 
side leng For some other approaches 


Gomerg |2003; [Hajian & Souradeep |2003; [Opher 
Dineen et_alj|2005; |Gomero et_al||2000; [Reboucas 
2000; 

In this work we analyze the most recent seven-year 
WMAP data for signatures of the three flat topologies 
of space mentioned above using the full covariance ma- 
trix of temperature-temperature fluctuations. By us- 
ing the symmetry groups of the spaces we construct 
efficient algorithms for the theoretical computation of 
the covariance matrix and the likelihood function us- 
ing that matrix. These algorithms can be used again as 
soon as the high precision data from the Planck satellite 
are released. The compu- 
tation of the covariance matrix is done using a modified 
version of the CAMB program (Lewis et atl 200d), as 
discussed in Sec. [B] and that of x? and the likelihood 
using the available WMAP code (Jarosik et al] 
. We have used 
only the TT correlations in our analysis. Including TE, 
EE, and BB correlations is straightforward, but would 
quadruple the computer time needed, without much im- 
provement in the results since these other correlations 
have much larger errors. 

There have been speculations in the literature 
about the detection of a special direction in the 
CMB map in which the first few multipoles of 
temperature-temperature correlations seem to be 
aligned (Land & Magueijo 
This is referred to as the “axis of evil” and is 
given by b = 60°, 1 = —100° in galactic coordinates. 
The topologies that we consider are not rotationally 
invariant, in particular T? x R!, and S1 x R? have one 
special direction (the infinite one in T? x Rt, and the 
finite one in S1 x R?), so we analyze the case where 
this special direction coincides with the axis of evil, to 
see if the axis of evil can be explained by one of these 


topologies. The authors of (Cresswell et _alJ{2006) have 


analyzed the topology S! x R? with the conclusion that 
it is not the explanation for the multipole alignment. 

The work is organized as follows. In section [| we 
present a slight generalization of Linde’s argument for 
the quantum creation of compact universes. We describe 
the calculation of the covariance matrix and the likeli- 
hood in sections [3] and [4] respectively. We present our 
numerical results in section [5] and discuss the goodness 
of fit, and maximum likelihood confidence intervals. We 
have done several checks of our analysis, which are given 
in Sec. [6] The possibility that our results are generated 
by a random fluctuation are analyzed in section[7] where 
we discuss Monte-Carlo skies. The possibility of spuri- 
ous effects due to a small residual CMB dipole in the 
data is investigated in Sec. [8] We summarize in section 
[9] Unless otherwise stated, everywhere in this work the 
side length of the torus is given in units of the distance 
to the last scattering surface Lo. 


2. QUANTUM CREATION OF COMPACT 
UNIVERSES 


Consider the standard Einstein-Hilbert action of grav- 
ity minimally coupled to matter. Here we will be only 
interested in “quantizing” gravity, so for the matter por- 
tion we will just consider energy density V without wor- 
rying about where it comes from. Then the action takes 
the form (ñ = c = 1, Mp, = (81G)~1/? = 1) 


1 
S= 5 | devar- 2V). (1) 
Following the standard procedure to derive the Wheeler- 
DeWitt equation we use the ADM form of the spacetime 
metric (Arnowitt et _al.|(2008) 
ds? = —N7dt? + hij(dx’ + N‘dt)(da? + Nidt). (2) 
The action can be rewritten in the form 


=de, (3) 


with 
L= “ (sR + a3 (By B — E*)- anv) (4) 
where °R is the 3-curvature of spatial slices, 
E;j = f = ViN; — V} Ni), 
E=£'. (5) 


There are numerous possibilities for the spacetime 
manifold and it may be described by infinitely many pa- 
rameters, so to be able to proceed we consider mani- 
folds with finite homogeneous spatial slices which can be 
characterized by one length scale a(t). In other words, 
we assume that locally the manifold is characterized by 
a Friedmann-Robertson-Walker metric while globally it 
can have any finite topology that is compatible with the 
metric. So by a suitable choice of coordinates we get in 
this case 


N=1, N,=0, (6) 
hij = a’ (t)kij ; (7) 


where the tensor kij is constant (it only depends on the 
choice of the manifold but does not depend on any of the 


U(a) 


a 


Fic. 1.— Plot of the potential U(a) for a positively curved space 
(8 > 0). The axes are in arbitrary units. 


coordinates). Then 


1\ 2 
Ej E“ — E? = —6 6) À (8) 
a 
Since we assumed a homogeneous spatial submanifold 
characterized by single length scale a, by dimensional 
analysis the volume must be proportional to a? and the 
curvature to a~?. Namely, 


[ee Vh=aa’, (9) 
3SR= > (10) 


where a and ĝ are dimensionless constants that depend 
only on the choice of the manifold. The Lagrangian then 
takes the form 


L = Š (aß — 6aa? — 2a°V) . (11) 


Now we treat a as the dynamical variable describing 


the geometry. The canonical momentum is then 


OL : 
Pa = a = —6aaa , (12) 


and the Hamiltonian becomes 


H=p,a—-—- L = —— (—p? — 6a? Ba? + 1207a*V) : 
12aa (13) 


Finally, we canonically quantize, replacing pa by the 
operator —i(d/da) to get for the Hamiltonian 


1 d? 
H = — | — = ba? Ba? + 12a?atV | x 14 
Daa (= 6a“ Ba" + 12a*a*V (14) 


Consider now the quantum creation of the universe 
with zero energy. Then the wavefunction of the universe 
W(a) satisfies the analog of the Schrödinger equation 
with Hamiltonian Eq. (14), which is called the Wheeler- 
DeWitt equation. In this case it takes the form 


d2 
— — 6a*Ba? + 1207a*V ) U(a)=0. (15) 
da? 
The effective potential energy is 

U(a) = 6a? Ba? — 12070*V, (16) 


and is shown in Fig. [I] which decreases to —oo for large 


4 


a since the second term dominates. However, for small a 
the first term dominates, so for 6 > 0 there is a po- 
tential barrier from a = 0 to \/G/2V, ie. the uni- 
verse has to first undergo tunneling before the expan- 
sion can start. This is the reason why the probability 
of quantum creation of positively curved spaces, which 
have 6 > 0, is thought to be highly suppressed com- 
pared to flat and negatively curved spaces. The action 
for tunneling through the barrier is S = aVV /3, where 
ao = 1/8/(2V) is the size of the created universe. The 
tunneling probability is x exp(—S), and is greater for 
smaller universes; S — 0,a9 > 0 as V > oo. For a flat 
universe, 3 = 0, and the barrier vanishes. 


3. COVARIANCE MATRIX CALCULATION 


Now we turn to the calculation of correlations between 
CMB temperature anisotropies in the flat topologies T°, 
T? x R!, and S! x R?. Locally they all look exactly 
like the infinite flat R? so Einstein’s equations and there- 
fore the Friedmann equations are unchanged from the 
infinite case. The calculation for the infinite case is de- 
scribed in standard textbooks (for a detailed derivation 
see (Dodelsonl [2003)), so let us briefly summarize that 
calculation and then focus on the differences between 
the infinite and finite universes. Essentially, one has to 
take the Einstein’s equations that describe the interac- 
tions between gravity and all of matter and Boltzmann’s 
equations for interactions between various types of mat- 
ter (most importantly, electrons and photons) and solve 
for the distribution of photons today given initial condi- 
tions set by inflation. Since the temperature anisotropies 
in the CMB are about five orders of magnitude smaller 
than the background, the calculation is done using per- 
turbation theory around the homogeneous background 
and keeping only first order terms. Then all of the differ- 
ential equations become linear and can be treated easily 
in Fourier space. This is where there is a key difference: 
in an infinite universe the spectrum of the Fourier modes 
k is continuous, while for compactified dimensions the 
spectrum becomes discrete. For a torus with side lengths 
In, Lə, and Ls we have k = (ki, k2, k3), 


2 2 
h= L m, hesi na ka= ns (17) 
L3 


Lı 
where nı, n2, ng are integers (the torus is essentially 
a box with periodic boundary conditions). So all of the 
equations in Fourier space remain unchanged, all we have 
to worry about is integrations over k which have to be 
replaced by sums 


d?k 1 
o ? Lı L2L3 pa (8) 


k 


over the discrete k values in Eq. (7). The set of points 
Eq. (£7) will be referred to as the k grid. 

The three cases studied here can be characterized by 
different values for L;. The three-torus T? has Ly = Ly = 
Ls = L, T? x R! has Ty = Lo = L, Ls =o, and finally 
St x R? has Ly = Lə = œo, L3 = L. All three cases 
can be treated in a unified manner by using the integral 


n SE 


dk bia 
Mimi'm’ = | gpro ja Yim (a 


notation, with the understanding that the integral is to 
be replaced by a summation if the corresponding L; is 
finite. 

The first set of summations over k arises when con- 
structing collision terms in Boltzmann’s equations. How- 
ever, we will not worry about these integrals for the fol- 
lowing reason. The Boltzmann’s equations are important 
only before the decoupling epoch, which corresponds to 
a redshift of about z ~ 1100. The comoving horizon at 
that time was about 50 times bigger than currently, and 
the current bounds on the size of the torus are of the or- 
der of the size of horizon, so at the epoch of decoupling, 
the size of the torus was at least about 50 times bigger 
than the causally connected part. As we will see later in 
section [5] the sums rapidly converge to the correspond- 
ing integrals when the topology scale is around 3 times 
the radius of horizon, which implies that the effects of 
finiteness can be safely ignored for the epoch of decou- 
pling (and before), and k can be treated as a continuous 
variable. All the equations are solved in Fourier space, 
exactly as for the infinite case. 

There is a summation over k when the final answer 
for the temperature fluctuations has to be converted 
from Fourier space back to real space. So let us pick 
up from that point in the calculation. The tempera- 
ture anisotropies O(n, x) in direction n at a given point 
x (chosen to be our location x9) are decomposed into 
spherical harmonics 


o0 l 
O(n, x) = 5 5 Qim(X)Yim (fi) ; (19) 


l=1 m=-l 


where the position space aj, coefficients are given in 
terms of the Fourier space temperature fluctuations 
O(f, k) by 


3 
dim (x) = | pet fan = (â)O(ñ, k). (20) 


The observed CMB fluctuations are given by the corre- 
lations between the different aj,,’s, 


Mimim' = (Aim (X0) Qj? m’ (xo)) z (21) 

The correlations between temperature anisotropies in 

k-space are related to the initial matter power spectrum 
(O(k, â)O*(k', n’)) 

O(k,k-f) O*(k,k-n’) 
ô(k) ô* (k) 
where the matter power spectrum is defined by 

(Dodelson 
(6(k)6*(k’)) = (27)?53 (k — k’) P (k) . (23) 


The ratios 6/6 on the right hand side of Eq. do 
not depend on the initial conditions since the equations 
are linear. All of the information about initial conditions 
is now absorbed into P(k). From Eq. (20), (21), and 
we get 


= (27)?6°(k — k’) P(k) , (22) 


OR Sat Yow EE, (24) 


ô* (k) 


Expanding O(k,k- ñ) into Legendre polynomials 
O(k, k- â) = X (~) (21 + DP(k-A)O(k), — (25) 
l 
and using the identity 
J dQ Py (Ê - â)Yim (Â) = 
we finally get 
Mimin = (4T)? (~i) x 
d3 *(k 
[ hrm eiw 
(27)° d(k) d*(k) 
We first review the standard result for an infinite uni- 


verse. In this case, the angular part of the integral over 
k in Eq. (27) can be done analytically giving 


Mimin = ô Ômm' Cı , (28) 


= Z / dk k? P(k) aoa 


The derivation remains the same in the finite case ex- 
cept that the k integral must be replaced by the sum 
Eq. (£8), so instead of Eq. (27) we get 


ÅT 4 
Yim (È), 2 
ea. îm (k) (26) 


Yia (ik) Yirm‘(k) (27) 


with 
2 


(29) 


Mimin = (4r)? (— i)li” x 
_ Sil) SEW yf 


Now we have to compute a three-dimensional sum 
Eq. (80) instead of a one-dimensional integral Eq. (29) 
which requires much more computational time. Also, we 
have to calculate all matrix elements with different l, m, 
l’, m’ whereas in the infinite case all 1 4 l or m 4 m 
(non-diagonal) elements vanish, while the diagonal ones 
do not depend on m. The reason for this is clear. In the 
infinite case, the problem has full rotational invariance, 
so that angular momentum is conserved. In the cases we 
consider, rotational invariance is broken. Even though 
rotational invariance is broken, there is still a large resid- 
ual discrete symmetry group which can be used to sim- 
plify the problem, and reduce the computational time. 
We will refer to this residual symmetry group as G. For 
the T? x R! and St x R? cases, G is the symmetry group 
of a rectangular parallelepiped with two sides equal, the 
tetragonal group D4 with 16 elements, whereas for T°, 
G is the symmetry group of the cube, the octahedral 
group Opn with 24 elements. 

The angular part in the sum in Eq. (80) can 
i separated (this has been suggested earlier in 


Phillips & Kogut]|2006)) 


4r)? (~ii oš 
Min m= SPO so a - 


ô(k 
x YO Yn (É)Yrm (k), (31) 


|k|=k 


where the first sum is over all the allowed spheres in the 
k grid while the angular sum is over a fixed sphere and 
depends only on the choice of that sphere. 

We can simplify the computation using the discrete 
symmetry group G of the manifolds 0,1,2. Consider a 


5 


fixed sphere of radius k. If the point (0, ¢) of the sphere 
is on the grid, then so are (6,6 + 7/2), (0, +T), and 
(0, ġ + 37/2). The angular sum over these four points is 
proportional to 


ém m) (1 4 eilm’—m) 5 + ilm —m)r 4 en) 


which is 0 unless m’ — m is divisible by 4, in which case 
it becomes 4e(""—-™, Consider the points (0,¢) and 
(—0,¢), which both lie on the sphere. Since 


Yim (—6, 6) = (-1)!- Vim (0, $), 


the sum over those two points is 0 unless 1 +l — m — m’ 
is even, but m + m’ is even if m’ — m is divisible by 
4, so the extra condition we get is that l’ —1 has to be 
even (this also follows from parity). The eight points 
+0,¢ + nT/2), n = 0,1,2,3 lie in the eight different 
octants, so the point (0,¢) can be chosen to lie in the 
first octant. To summarize, the angular sum is nonzero 
only if l’ — l is even and m’ — m is divisible by 4, in which 
case it is equal to 8 times the sum over one octant. Extra 
care is needed for points on the boundary of the octant 
to avoid double counting. 

Consider the points (0, ¢) and (0, 7/2 — ) correspond- 
ing to swapping nı with ng. Taking into account that 
m’ — m is divisible by 4, we get 


eilm’—m)¢ a im =m) (5-4) = 2 cos ((m! _ m) $) , 


which implies that the angular sums are real. Further- 
more, (—i)!i" is also real for even I’ — l and P(k) and 
1(k) /ô(k) are real, so the covariance matrix elements 
Mimvm are all real implying Minium: = Mimim. Also, 
since Y; (0, ¢) = (-1)'Y,*, (0, @) and m’ — m is divis- 
ible PY 4, we get Mimim = Mi, -m , -m . 

T? has more symmetries which can be used to further 
speed up the calculation for this case. For example, the 
sums in Eq. (31) over the spherical harmonics are the 
same for all L. Changing L is a rescalling of the allowed 
momenta by 1/Z. Thus the angular sum for |k| = 
for a T? of size L is the same as the angular sum for 
|k| = Ak for T? of size L/A. Thus the angular sums can 
be computed once, and then used for all values of L. 

Since the calculation of 0;(k)/d(k) is identical to the 
case of infinite flat universe, we use the well-known 
FAST (Seljak & Zaldarriagal ) for that part of the 
calculation. Our code is essentially a modification of 
CAMB. It takes the sides of the torus as extra input 
parameters and outputs not only C; but also the com- 
plete matrix Mimim:. For large values of l, the discrete 
sums over k approach the continuum result, so we only 
use Eq. for | < 30, and use the continuum result for 
l > 30. The difference between the discrete and contin- 
uum values for Mjmim is less than 0.5% for l = 30. As 
an example, in Fig. 2] we have plotted the ratio of the 
power spectrum C; for Mı = T? x R! with L/Lo = 1.9 
to that for infinite space R?, where C, has been defined 
as 


~ 


l 


1 


for Mı = T? x R!. The l = 2 power is reduced by 
7%, and the ratio of power spectra oscillates and rapidly 
approaches unity. The two differ by 0.1% at l = 30. 


C,/C@) 


0 10 20 30 
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Fic. 2.— Plot of the ratio of C, for 4 = T? x R! with L/Lo = 
1.9, to that for infinite space R3. 


4. LIKELIHOOD CALCULATION 


The matrix Mim'm’ computed as discussed above is 
compared to the experimental data from the 7-year 
WMAP survey. Since rotational invariance is broken, we 
need to vary the orientation of the torus relative to axes 
fixed in space, to find the best fit. We do this by rotating 
the data relative to the torus in computing the likelihood 
function. We specify the orientation by three Euler an- 
gles (¢,0,w) in the following way. The axes x,y,z are 
fixed in the coordinate frame of the CMB data, i.e. the 
observed universe, and the x’, y’, z’ axes are fixed in the 
torus. Start with the CMB-fixed and torus-fixed axes 
aligned. Rotate the torus counterclockwise around the 
z-axis by angle ¢, then around the new x-axis by angle 
0, then around the new z-axis by angle ~ to get the final 
torus orientation. The angles ¢ and 0 give the orienta- 
tion of the z-axis of the torus while the angle Y gives the 
orientation of the torus around its z-axis. We can make 
use of the symmetries of our topologies to speed up the 
calculation since various Euler angles can give equiva- 
lent orientations of the torus. Two sets of Euler angles 
(, 6, Y) and (¢’, 0’, w’) are equivalent if 3g € G such that 


Ro, 0,0) = R(g) RO, 4") (33) 


where R(¢,6,) is the coordinate transformation rota- 
tion matrix corresponding to (¢,@,wW) and R(g) is that 
corresponding to the discrete element g. This defines 
an equivalence relation on the set of all possible Euler 
angles. We take a uniform grid on all possible angles, 
then divide that grid into equivalence classes according 
to Eq. and take one representative from each class. 
We have scanned over ~ 4000 inequivalent angles. 
Different orientations of the torus were considered 
in the previous analysis of first-year WMAP data for 
T? (Phillips & Koguti[2006), but they only considered a 
uniform grid on the range 0 < ¢,0,w < 2/2. Note that 
this does not cover all possible orientations of T3. The 
first two angles describe the orientation of the z-axis and 
by their assumption on the range of ¢, 0, ~, the z-axis al- 
ways lies in the first octant. However taking into account 
all the symmetries of the cube there are 6 equivalent axes 
that can play the role of the z-axis, the +a, +y, and +z 
axes, while there are 8 octants. In other words, there are 


possible orientations of the cube for which none of the 6 
axes lies in the first octant. 

After choosing a torus orientation, we calculate the 
likelihood in the real space of orientations on the last 


scattering surface. For N, pixels the likelihood function 
is given by? (Dodelson 


1 


Pn 
(20) edet C2 P 


(-5aro=a) , (34) 


and the y? function by 
xX = ATCA, (35) 


where A; is the vector of pixels and Ci; is the covariance 
matrix that also includes the noise. The indices 7, 7 label 
the different pixels, which are in directions ñ; j on the 
sky. 

We now describe how to convert the matrix Mimim to 


the covariance matrix C;;. The temperature fluctuation 
measured in a pixel i is given by (Dodelson! 
6; = fan O)B:(à), (36) 


where B; is the beam pattern at the pixel į and is specific 
to the experiment. Usually the beam patterns have the 
same shape for every pixel and are axially symmetric 
around the center of the pixel, as is the case for WMAP, 
so if we denote the direction to the center of the pixel 
by ñ; then the beam pattern can be decomposed into 
spherical harmonics 


Bi (å) = X BiYim (fi) Yin, (A) - (37) 


lm 


Using Eq. to decompose ©(n) into spherical har- 
monics, we get for the theoretical covariance matrix 


Ci; = (0:9;) = 5 Mimi m! BiBr Yim (ti) Yim (yj) - 
lml’m!’ 

(38) 
In computing C;;, we have to vary the orientation of 
the torus relative to the sky. In implementing the Euler 
angle rotation, one can compute the Mjmym matrix in 
the torus-fixed coordinate system, so that it remains un- 
changed as the Euler angles are varied. The pixel direc- 
tions n; are changed to n; > R(¢, 6,1) n;. Equivalently, 
one can work in the CMB-fixed coordinate system, and 
rotate the torus, which gives Mjmym transformed by the 
angular momentum rotation matrices, 


Mimi'm! => 5 Mmrw DË*(R)D® (R) š (39) 


n 
n,n! 


Note that in the infinite universe case, Eq. holds, 
and the result Eq. simplifies to 


At eR 
Ci; = a Bi O Piu Ày), (40) 
l 


independent of the rotation R(¢, 0, w). 

The computation of the covariance matrix using 
Eq. is more involved than the infinite case, Eq. (40), 
so the likelihood calculations require far more computer 


2 To distinguish from length scale L, the likelihood is denoted 
by @ everywhere in this work. 


time than the conventional case. C;; must be recalcu- 
lated for each set of Euler angles. There are 458403 in- 
dependent elements in Mjmym for 2 < l < 30 of which 
57840 satisfy the | = l’ (mod2),m = m’ (mod 4) condi- 
tion, and 2482 values for each of the indices ¿i and j. The 
slowest step in the computation is evaluating the sums 
on l,m, l,m’ in Eq. (88) for all values of {7,7}. 

An Euler angle rotation of the sky maps points on the 
sphere to rotated points on the sphere. For infinitesimal 
pixels, this corresponds to a reshuffling of the pixels, i.e. 
if pixel 7 at n; is mapped by the rotation to nj, then pixel 
i — pixel j. An exact reshuffling of pixels would greatly 
simplify the computation — instead of recomputing C;;, 
one could simply permute the indices on Ci; to get the 
transformed matrix. In particular, det C would remain 
invariant under this transformation. 

The WMAP pixels have been chosen using the 
HEALPix grid (Gorski et al]/D005). The pixels are cho- 
sen to lie along lines of constant lattitude, and they have 
equal solid angles. This implies that the spacing of the 
pixels varies as a function of lattitude. As a result one 
cannot treat rotations of the sky as a pixel reshuffling 
transformation. One can approximate the rotations by 
a pixel transformation by mapping the rotated pixel to 
the one closest to it in the HEALPix grid. The likeli- 
hood computed using this method differs from the exact 
result using Eq. (88), and is not accurate enough for our 
purposes. The above approximate relation between rota- 
tions and pixel permutations does, however, explain why 
det C is approximately independent of the Euler angles. 

For a finite universe, one has to use Eq. (88) with 
the value for Mijmym: computed as described in Sec. [B] 
The finiteness of the universe only affects the large-scale 
anisotropies, so the difference between the infinite and 
finite cases goes to zero with increasing l. For that rea- 
son we will look only at low-l portion of anisotropies, 
l < 30, and use the infinite manifold result Eq. (40) 
for | > 30. We calculate y? and the likelihood -Z us- 
ing a modification_of the likelihood code provided _b; 
the WMAP team (Jarosik at al 2011 Larson et al BOL 
as a function of the new parameters 
(L, ¢, 0, Y). Since we are interested only in low-/ effects 
we use the low-resolution portion of the likelihood code. 
We use the experimental data in the exact same form as 
provided by the WMAP team without any further mod- 
ifications. The temperature map used is the smoothed 
and degraded ILC map with the Kp2 mask applied to 
remove the galactic plane and strong point sources. The 
map originally has 3072 pixels, but only 2482 are left 
after the mask. 

Ideally, one would have to do a fit to the experimen- 
tal data varying the four new parameters (L, ¢, 0, y) in 
addition to all the other cosmological parameters. The 
cosmological parameters affect the whole spectrum of 
anisotropies while only the low-/ part of the spectrum is 
affected by the new parameters, so we fix the other cos- 
mological parameters at their best-fit values as given by 
Misi soven gear NI Gale and 
only vary the new parameters. The values of the cos- 
mological parameters that we use are 
100Q,h? = 2.227, Qh? = 0.1116, Q4 = 0.729, 
Ns = 0.966, r = 0.085, A3 (0.002 Mpc~') = 2.42 x 107°. 
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Fic. 3.— Plot of Indet C/C; against the Euler angle ~ for fixed 
o, 9 for the .@ topology at L/Lo = 1.8, 


1073 


je) 


5. RESULTS 


We have computed the likelihood and x? for the three 
cases, so = T’, ™ = T? x R! and Mo = S! x R?, 
for different values of L/Lo as a function of the Euler 
angles. L/Lo ranges between a minimum value of 0.5 — 
1.0 depending on the manifold, and a maximum value of 
L/Lo = 2.6, in steps of L/Lo = 0.1. By the time L/Lo = 
2.6, the results are almost identical to the flat-space case 
Mo = RÌ. For the statistical analysis discussed later in 
this section, we have used interpolation to construct a 
smooth function of L. 

The relation between likelihood and y? is 


-2m Y= ? + Indet C/C; + Indet(2xCy) (41) 


where Cy is a fiducial covariance matrix used by the 
WMAP collaboration. Cp is independent of L and the 
Euler angles, and drops out of all likelihood ratios. x? 
and —2In& differ by Indet C/Cy (up to an irrelevant 
constant). The likelihood for s, three-dimensional flat 
space, will be denoted by -%., and is YY = 3573.4 

As noted earlier, for fixed L/Lo, IndetC/Cy varies 
weakly with the Euler angles. In Fig. [3| we have plotted 
the variation of ln det C/Cy as a function of 7, for fixed 
values of the ¢,@, at L/Lo = 1.8. The overall variation 
of Indet C/C} against W is less than unity. 

—2In¥ (and hence x?) has a strong variation with Eu- 
ler angles at fixed L/Lo. In Fig. [4] we have shown plots 
of the variation of 2In-% — 2ln Z with Euler angle w 
for fixed ¢,@ for the .4 topology. æ is independent 
of the Euler angles. The solid red curve has been chosen 
to have L/Lo = 1.8, and ¢,6 values that maximize the 
likelihood at this value of L/Lo. There is a large varia- 
tion of —2ln Z with the remaining Euler angle Y, and 
the global minum of —2 ln X is 2In Zs —2ln Z = —-17.2 
at w/(2m) ~ 0.05. The strong dependence of —2 ln X on 
orientation makes it difficult to find the true global min- 
imum of the —2In-¥ and x? functions. We have done a 
scan over all Euler angles with a spacing of 0.057, to iden- 
tify valleys, followed by a finer scan to find the mininum. 
By comparing our numerical minimum with the next best 
point, we can estimate the uncertainty in our minimum 
—2InY and x? values at less than 0.5. The dashed blue 
curve in Fig. [4Jis also for L/Lo = 1.8, but with ¢, 0 fixed 
at random values, rather than those for which —2 In Z 
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Fic. 4.— Plot of 2In Yoo — 2In Y against the Euler angle ~ for 
fixed ¢,6 for the o topology. The solid red and dashed blue 
curves are for two different values of 4,6 at L/Lo = 1.8, and the 
dotted green curve is for L/Lo = 2.2. The solid red curve is for 
L/Lo = 1.8 with ¢,6 fixed to be the best fit values, the dashed 
blue curve is for L/Lo = 1.8 with ¢, 0 fixed in a random direction, 
and the dotted green curve is for L/Lo = 2.2 with ¢, 0 fixed to be 
the best fit values. 


vs. w passes through the global minimum. There is still 
considerable dependence as one varies the third angle Y, 
but the dependence is much weaker than for the solid red 
curve. The dependence of —2 lIn Z drops rapidly with in- 
creasing L/Lo. For L/Lo = 2.2, the dotted green curve 
in the figure, the overall variation is about 6.5. 

The plot of x? and —2InY against L/Lo is given in 
Figs. [5] [6] F] for the three cases, o = T?, 1 = T? x R! 
and % = S' x R?, respectively. We have plotted the 
maximum and minimum of x? over all possible orienta- 
tions of the torus at each value of L. There is a sig- 
nificant variation in y? as a function of orientation, as 
noted earlier, and some orientations are strongly pre- 
ferred over others. In each plot, x? ranges between 
the uppermost and lowermost solid black curves, as one 
varies the orientation of the manifold by varying the Eu- 
ler angles (¢, 0, Y). For the smallest values of L/ Lo, Ax? 
between the worst and best orientations is 377, 191, and 
156 for 4@,1,2, respectively. As L/Lo increases, the ef- 
fect of a compactified direction decreases. By the time 
L/Lo = 2.6, the fit results are very close to the case of 
the infinite manifold R°, and Ay? < 4 for the different 
orientations. 

We have been unable to find any pattern to the best- 
fit orientation ¢, 6, of the torus as a function of L/Lo. 
We have examined the possibility that the manifolds we 
consider are aligned along the axis of evil. To do this, we 
have chosen the preferred axis of the manifold (the z-axis 
for T3, the R direction for T? x R! and the St direction 
for St x RÊ?) to point along the axis-of-evil direction b = 
60°, l = —100° in galactic coordinates, and allowed for 
arbitrary rotations of the manifold around this direction. 
All the angles are varied with step 7/100 = 1.8°. This 
gives a subset of all the orientations we have considered, 
and the x? range has been plotted as the dashed colored 
curves in the figure. The colored curves lie between the 
black curves (as they must), but they do not lie towards 
the best-fit x? line. This shows that there is nothing 
in our computation that picks out the axis-of-evil as a 
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L/L, 
Fic. 5.— Plot of x? and —2In¥ against L/Lo for Mo = T? 


for different Euler angles (4, 8, v) The lower solid black curve 
(solid triangles) is the minimum x? (or —2In.#), and the upper 
solid black curve (open squares) is the maximum x? (or —2 ln 2). 
The lower dashed colored curve (solid triangles) and upper dashed 
colored curve (open squares) are the minimum and maximum y? 
(or —2In@) with a symmetry axis of the manifold restricted to 
point along the axis of evil. 


C.L. Mo Ma Mz 
68% 1.47 1.24 0.92 
90% 1.33 1.02 0.61 
95% 1.24 0.96 0.56 


TABLE 1 
LIMITS ON L/Lo USING THE x? GOODNESS-OF-FIT TEST. VALUES 
LESS THAN THOSE IN THE TABLE ARE EXCLUDED AT THE 
CONFIDENCE LEVEL GIVEN IN THE FIRST COLUMN. 


preferred direction. 

We use a goodness-of-fit test to see whether we can 
rule out the hypothesis that the universe has topol- 
ogy M; of size L. The minimum y? values are x? = 
2469, 2467, 2472 at L/Lo = 2.1, 2.1, 2.2 for Aoi, 
respectively. Previous studies have aso pound ae eee found an indica- 
tion of a dip in y? around L/Lo ~ 2.1 ( 
Using the computed values LEEI x? we have the 
EN given in Table [I] 

We estimate confidence intervals for L/Lo us- 

likelihood ratios Wash 
97D) (Nakamuract all DOU $33). Maximum ikeli 
hood confidence intervals are exact if the distribution 
is Gaussian. For non-Gaussian distributions, they have 
an error of order 1/,/n, where n is the number of data 
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Fic. 6.— Plot of x? and —2ln ¥ against L/Lo for Mı = T? x R! 
for different Euler angles (¢, 0, Y). See the Fig. B]caption. 


points [1938). The maximum likelihood value 
is at L/Lo = 1.8,1.9,1.9 for 41,2, respectively. The 
confidence intervals for L/Lo are determined by using 
Aln Y, the difference of ln- from the value that max- 
imizes the likelihood function 
(1971). Plots of the L/Lo confidence intervals as a func- 
tion of 1 — a, where a is the confidence level, are shown 
in Figs. [8] Q| 0] The data shows a preference for a finite 
universe with size L/Lo ~ 1.9 corresponding to L ~ 27 
Gpc. The allowed L range extends to L/Lo > 2.6 at 
a confidence level a = 1074 for Mlo, 2 x 107° for M 
and 4 x 1073 for .@. Thus the data show evidence for 
a finite universe at a confidence level a = 2 x 107° for 
the T? x R! topology. The 95% confidence intervals 
are L/Lo € [{1.7,2.1],[1.8,2.0],[1.2,2.1] for 0,1,2, 
respectively. 

We have scanned over ~ 4000 different orientations for 
each value of L, with a finer scan near the minima, so 
that the error in y? and —21In-Y is < 0.5. The difference 
in —2 ln Y between its value at L > oo and its minimum 
value (which occurs at L/Lo = 1.9 for .@) is 20.4, which 
is well outside possible numerical errors. Note that the 
main numerical uncertainty is finding the true minimum 
of —2ln# for finite values of L. The minimum value of 
—2ln¥ has been determined with an accuracy < 0.5. 
The actual difference in likelihoods between finite and 
infinite L can only be greater than what we have found. 
There is an indication that a finite universe fits the data 
better than an infinite one. However, the “standard” 50- 
criterion for a discovery, corresponding to a confidence 


DRT rey rr 
pre 


= 
1 


3500 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 
0.5 1 1.5 2 2.5 


aie 


Fic. 7.— Plot of x? and —2In ¥ against L/Lo for M2 = St x RÊ? 
for different Euler angles (¢,6,7). The two colored curves are 
almost on top of each other. See the Fig. [5] caption. 
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Fic. 8.— Plot of the L/Lo confidence interval as a function of 
the confidence level for T’. 
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level a = 5.7 x 1077, includes the value L = oo. 
The Euler angles for the best fit case 4 with L/Lo = 


1.9 are (@ = 21° 4 
which corresponds, 


1 2°. 0 = 53° 4 


E 2° = 61° + 2° 


for the infinite direction, to (b = 


37° 4 


-2° ] = 291° + 


2°) in galactic coordinates and (a = 


182°4 


2°, § = —25°4 


t 2°) in J2000 equatorial coordinates. 


This is close to the direction (b = 30° 42°, 1 = 276° 43°) 
of the velocity of the Local Group inferred from the CMB 
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Fic. 9.— Plot of the L/Lo confidence interval as a function of 
the confidence level for TŽ x R?!. 
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Fic. 10.— Plot of the L/Lo confidence interval as a function of 
the confidence level for St x R?. 


dipole (Lineweaver]1996). We discuss the possibility that 
our signal is due to a dipole contamination in Sec. [8] 

We can compare the probabilities that the universe 
has topology -4,1,2,3. We use the Bayesian prior that 
there are four discrete choices, @% with L/Lo = 1.8, Mi 
with L/Lo = 1.9, @ with L/Lo = 1.9, or the infinite 
case læ, each of which is equally probably. Since we 
are comparing four discrete cases, there is no ambiguity 
due to choice of measure in choosing an equiprobable 
prior. Then using the likelihood ratios gives the posterior 
probabilities 


p(M) =0.96, 
p (Me) =2 x 10-5, 


p(Maæ)=1.3 x 107°. 


The probability of the infinite universe Mls is very small, 
and the most probable topology is %4 = T? x R', with 
two dimensions compactified, and one infinite. 


5.1. Fisher Information 


The Fisher information can be used to compute the 
variance V of the length L determined using the maxi- 
mumn likelihood method. The Fisher information is given 


by 
a Ay (ony aa 0 
mesna e TL 


ƏL ƏL’ 
Using the covariance matrix C for the 41 toplogy with 
L/Lo = 1.9 gives 


(42) 


V-t =3.3 x 10° (43) 


so that the error estimate for L/Lo is VV = 0.017. The 
Fisher information error Eq. (43) corresponds to using a 
quadratic approximation to the likelihood function about 
its minimum to determine the error, and gives a smaller 
error than that obtained earlier using the exact likelihood 
function. 


6. CHECKS 


We have been unable to find a simple explanation for 
the better fit due to a finite topology. However, there are 
some possibilities which we can test. 

The measured cosmic microwave background 
anisotropy has a smaller value for the quadrupole 
power Cə than the theoretical expectation value. There 
is a large cosmic variance in C2, so this is not a discrep- 
ancy between theory and experiment. Fig. 2]shows that 
the predicted value of C2 for a finite universe is reduced 
from the infinite universe value. The greater likelihood 
for a finite universe is not due to lowering the value of 
C2. We have checked this by determining the likelihood 
using Mimim’ for the finite case, but with the l = l’ = 2 
values replaced by their values for the infinite universe. 
For Mı with L/Lo = 1.9, —2InY increases by 0.43, 
which is much less than the 20.4 difference in —2InY 
from the infinite universe. 

As another test, we have computed the likelihood for 
the “ topology for L/Lo = 1.9 by using a truncated 
Mimtm’ Matrix. The truncated matrix is constructed by 
using Mimym: for the finite topology for 5 < l,l < 20, 
and using Mimim’ for the infinite universe, i.e. Mimim: = 
Crow mm’, for | and or l’ outside this range. A plot of 
the likelihood as a function of the Euler angle ~ for this 
truncated Mjmim is plotted as the dashed blue curve in 
Fig. This can be compared with the likelihood curve 
using the full Mimim for the finite topology, shown as 
the solid red curve. The dip in the likelihood difference 
to —20.4 is the signal that the finite topology is a better 
fit than the infinite universe. The plot for the truncated 
matrix is similar to that for the full matrix, except that 
the small-angle fluctuations have been smoothed out, as 
is to be expected since higher / terms have been dropped. 
Note that the dip in 2In Z» — 2ln# is very similar in 
both cases, and the minimum of 2 In .%,,—2 ln -X is nearly 
the same. This shows that the effect we find is not due 
to the low-l modes (quadrupole, octupole), and is also 
not an edge effect as a result of only using l < 30 in 
the compuation. For 5 < l,l’ < 20, Mimim’ has 86736 
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Fic. 11.— Plot of 2ln “% — 2ln Y against the Euler angle w 
for fixed ¢,6 for the .@ topology for L/Lo = 1.9. The solid red 
curve uses the full matrix Mimi’m’, the dashed blue curve uses the 
matrix truncated to 5 < l,l’ < 20, and the dotted green curve uses 
the matrix Mimim’ Oy, retaining only the part diagonal in l. 


elements of which 11056 satisfy the | = l’ (mod2),m = 
m’ (mod 4) condition and are non-zero. 

The off-diagonal elements in Mjmym’ are important for 
the calculations. We have also computed the likelihood 
by retaining only the elements which are diagonal in J, i.e. 
using Mimms. This drops the elements in Mimym/ dw 
which are off-diagonal in l, while retaining the elements 
which are off-diagonal in m for a given l. The likelihood 
with this matrix is the dashed green curve in Fig. 
With this matrix, the likelihood deviates much less from 
the infinite universe, and the dip near ¢)/(27) ~ 0.17 is 
much less pronounced. 


7. MONTE-CARLO SKIES 


The results of the previous section were obtained us- 
ing a likelihood analysis of the WMAP7 data. One can 
study whether the better fit of a finite topology is due 
to a statistical fluctuation. Since the big-bang is not a 
repeatable experiment, this must be done by generating 
random Monte-Carlo data for the pixels A;, and redo- 
ing the analysis for this Monte-Carlo data. To actually 
do this numerically is beyond the computing power we 
have available. Luckily, for the problem at hand, we can 
analyze the Monte-Carlo problem analytically. 

Assume that the pixels A; are generated by the co- 
variance matrix Cə for an infinite universe, so that the 
probability distribution is 


ee ee 
PA) = Taane »( z^ ozaj- (44) 


The likelihood function computed using A; and covari- 
ance matrix C (of a finite universe) is 


— 2log -Z = ATOT!A + Indet(2aC), (45) 


and the likelihood constructed using the covariance ma- 
trix Cs of the infinite universe is 


— 2log -Z = ATOZ!A + lndet(2nC). (46) 
Let 
h=(—2log Z») — (—2log 2), (47) 
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be the difference of the two log-likelihoods. In our anal- 
ysis, we found h = 20.4 > 0, so that the finite universe 
was more likely than the infinite universe. The average 
value of h over Monte-Carlo data can be computed using 
Eq. and Eq. (47). The two-point function is 


(Aj;A;) = (Cæ) (48) 


ij ? 
so that 
lh) =N — Tr C!C + Indet(C..) — Indet(C) .(49) 
It is convenient to define the symmetric matrix 
S = C1420710, (50) 


which is a positive matrix since C and Cy are positive 
matrices, and has eigenvalues s; > 0. In terms of S, 


(h) =N — Tr S + Indet S 
=ġ_ [1 - s; +lns;] . (51) 


The function 1—s+Ins < 0 with its maximum at 0 when 
s=1. Thus 


(h) <0 (52) 


and (h) = 0 only if S = 1, ie. C = Cy. This gives the 
intuitively obvious result that the best fit for data gener- 
ated with covariance matrix Cə is, on average, given by 
fitting using the same covariance matrix Cs. Any other 
covariance matrix C used for fitting, on average, gives a 
lower likelihood. 

If instead of Eq. (47) we had used the difference of x7, 


hy =x% E x’, (53) 
then 
(h) =N = The "Ce =N-TrS=)5° [1-5] ; 
(54) 


and (h,) could have either sign, since s; > 0, but need 
not be smaller than 1. For example, a simple rescaling, 
C = àC», with à > œ can always make y? — 0, its 
minimum poissible value. This option is eliminated for 
likelihood because of the det(27C) term. 
Using for C the best-fit topology Mı with L/Lo = 1.9, 
we find the numerical values 
N =2482, 
In det (Cso) = 1097.8 + Indet Cy, 
In det(C’) = 1082.73 + Indet C} , 
Tr C!C = Tr S = 2516.6, (55) 


so that 
(h) =—19.5. (56) 


This differs from the value we find of h = +20.9 by Ah = 
h — (h) = 40.4. The probability that Ah is a statistical 
fluctuation can be determined by computing the variance 
of h using the four-point function 


(A; Aj Ax Az) = (Coo) ij (Coo) it T (Coo) ix (Coo) 51 
+ (Coo) 54 (Coo) jx > (57) 
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to obtain 
((Ah)?) =2N +2TrO-10,,0-1 C55 — 4 (Tr) 
=2Tr (1— C710.) =2Tr (1-S)?. (58) 

In our case, 

TG 'C..0 *C.—26100, (59) 
so that 

((An)’) = 117.6 = (10.8). (60) 


Our observed value of Ah = 40.4 is 3.70 away from 
the mean, so the probability that a fluctuation gives h 
larger than or equal to our observed value is 1.1 x 1074, 
assuming a normal distribution. 

While the distribution of the data A; is Gaussian, the 
distribution of the likelihood difference h is no longer 
Gaussian. We can also compute higher order connected 
correlation functions of h, 


((Ah)"),=2" (r — 1)! Tr (1-C10,,.)" 
=2" (r—1)! Tr A-S), (61) 


from the generating function 
log (e*4) = —XTr (1 — C7* Coo) 
1 
= it ln [1 — 2A(1 — C71 Ca )] 


=-XIi (1-8) -5Tr In [il — 2A(1— 5)] , 


(62) 
so that 
((Ah)?) =8 Tr (1 - 071C)? 
((Ah)*) =48 Tr (1-C710,,)* , (63) 
where the fourth-order correlation is 
((Ah)*) = ((Ah)*), +3 (Ah)? , (64) 


in terms of the connected correlation. The mean value 
(h), and all the connected correlation functions ((Ah)")., 
are of order N, the number of data points. Thus the 
relative correlation ((Ah)") / (h)" is of order NH=". 

The numerical values for our case are 


((An)°) = 489.0, 
(aw) = 4303.7. (65) 


We can get a better estimate of the probability that 
h = 20.4 is due to a statistical fluctuation by using these 
higher order moments. We have fit (h) and ((Ah)”), 
r = 2,3,4 to a probability distribution 


p(h) = po exp[—(h — ho)? — c2(h — ho)? 
c3(h — ho)? — ca(h — ho)*] , (66) 


and found using this distribution that the probability 
that h — (h) > 40.4 is 1078, which is smaller than the 
value obtained earlier using a normal distribution for h. 


8. DIPOLE CONTAMINATION 


The symmetry axis of .4, points in the direction 
(b = 37° 2°, 1 = 291° + 2°), which is close to the direc- 
tion of the velocity of the Local Group (b = 30° + 2°, l = 
276° + 3°) (Lineweaver [1996). The CMB has a large 
dipole asymmetry of 3.358+0.001+0.023 mK in the direc- 
tion (b = 48.05° + 0.11°, 1 = 264.31° + 0.2°) (Lineweaver 
[1996). Suppose that the data is contaminated by a dipole 
contribution that has not been properly subtracted out.? 
Could a residual dipole explain the results we have 
found? 

To study the effect of a residual dipole, assume that 
the observed pixels are 


Ass = A; +p-f; = A; + di, 


d;=p-fi; (67) 
where A; are the true fluctuations given by the distribu- 
tion Eq. and p is the residual dipole contamination 
in the data. Then Eq. (45146) are replaced by 
— 2log Y=(A +d)" C7! (A + d) + Indet(2zC) , 
—2log Z, = (A + d)? Cr! (A +d) + ln det (27C) . 
(68) 


From these, we find 


(h) =N — Tr C!C + mdet(C,,) — ln det (C) 
Co d) (a Ca), 
(Ah?) =2N — 8d?C7'd + 4d" Otd + 4d" CCC 7d 
+2Tr C7'Cx.C~'Cx, — 4(TrC7*C.) . (69) 
Dipole contamination produces a systematic shift in h 


from its value in Eq. given by the (d?’CZ/d) — 
(dT Ctd) terms, which can be written as 


(d? Cx 'd) — (dOd) = paps Dag (70) 


in term of the components Pa = (Px, Py, Pz) of the dipole. 
We find 


1.22 —0.004 —0.114 
D = | —0.004 —0.0281 —0.003 
0.114 —0.003 —0.0127 


mK? (71) 


The largest eigenvalue is —1.23mK~?. To get a shift in 
—2lnY of 20.4 requires a dipole contamination |p| of 
around 4mK~?. This is larger than the observed dipole, 
and several hundred times the quoted uncertainty in the 


CMB dipole (Lineweaver|1996), and is excluded. 


9. CONCLUSIONS 


We have analyzed the possibility that the universe 
has compact topologies % = T3, æ = T? x R! and 
Ma = S' xR? using modifications of the available CAMB 
and WMAP 7-year likelihood codes. The maximum like- 
lihood 95% confidence intervals are 1.7 < L/Lo < 2.1, 
1.8 < L/Lo < 2.0, 1.2 < L/Lo < 2.1 for the three cases, 
respectively. Using the Bayesian analysis discussed ear- 
lier, we find that the most probable universe has the 
compact topology .4. An infinite universe is compati- 
ble with the data at a confidence level of 2 x 1075 (i.e. 
4.30). We find no evidence of a preference for the axis- 
of-evil direction. The improved fit for a finite universe 


3 This possibility was suggested to us by B. Keating. 


is not due to the lowered prediction for the quadrupole 
anisotropy; this accounts for only a small fraction of the 
increase in likelihood. 

It would be useful to investigate whether any system- 
atic effects in the WMAP data introduce effects with 
cubic symmetry that can mimic the effects of a torus 
topology. The best fit results do not pick out any special 
orientation for the torus, such as a torus with symme- 
try axis perpendicular to the galactic plane, that might 
lead to systematic effects that lead to a fake signal. Pix- 
elization of the data using the HEALPix grid also should 
not introduce cubic symmetry terms along an axis not 
aligned with the galactic pole. The best fit results for 


13 


Mı have symmetry axis which is near (~ 10°) the direc- 
tion of the Local Group velocity. 
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